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Abstract 

A system of N bosons in a two-dimensional harmonic trap is considered. The system 
is treated in term of the slave boson representation for hard-core bosons which is 
valid in the arbitrary density regimes. I discuss the consequences of higher order 
interactions on the density profiles by mapping the slave boson equation to the 
known Kohn-Sham type equation within the Density functional scheme. 
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1 Introduction 

The possibility of Bose-Einstein condensation (BEC) phase transition occuring 
in a 2D geometry has a remarkable history and breakthroughs. Recent exper- 
iments have realized a quasi-two-dimensional (Q2D) cold Bose gas by tuning 
the anisotropy of the trapping potential [1,2]. This achievement has attracted 
vast attention on the new physics appearing in the currently studied low- 
dimensional trapped gases and the nature of finite-size and finite-temperature 
effects. The effect of dimensionality on the existence and character of BEC 
or superfluid phase transition has seen a spurt of interest in recent years. A 
peculiar feature of the low-dimensionality (ID or 2D geometry) at finite tem- 
perature is the absence of a 'true condensate' following the Bogoliubov k~ 2 
theory, believed to be originated from the long wave-length phase fluctuations 
[3] . This has been proved in the seminal work of Mermin and Wagner [4] and 
Hohenberg [5]. A Kosterlitz-Thouiess transition to a superfluid state still oc- 
curs on cooling through the binding of vortices of opposite vorticity [6] and 
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leads to an algebraic decay of the one-body density matrix. However, in a 
trapped 2D fluid the modification of the density of states caused by the con- 
fining potential allow a true condensate to exist even at finite temperature. 
Another important consequence lower dimensionality is that the T-matrix for 
two-body collisions in vacuo at low momenta and energy in the dilute regime 
vanishes [3,7]. It is then necessary to evaluate the scattering processes between 
pairs of Bose particles by taking into account the presence of surrounding con- 
densate particles [8-11]. 



The dilute limit (pa 2 <C 1), where p is the uniform density of particles and a is 
the s-wave scattering length of a Bose gas can be described quite satisfactorily 
by traditional approaches, such as the Gross-Pitaevskii(GP) approach in the 
weakly interacting regime [12]. However, it is intuitively obvious that as the 
density of particles in the condensate or the scattering length increases, the 
nature of interactions will become more and more complex, and one will need 
to take even three-body and higher order effects into account. It is a crucial 
point since in the strictly two-dimensional trap, the scattering length of the 
bosons become larger than the axial thickness of the condensate cloud. Thus 
we need to go beyond the mean-field approach to describe the system of dense 
Bose gas . 



These limitations motivates me to address the problem particularly the 2D 
dense condensate by introducing the slave boson model (SB) which is valid for 
arbitrary density regimes. The slave boson technique was first introduced by 
Kotliar and Ruckenstein [13], who used it to deal with the fermionic Hubbard 
model. A functional integral approach to the problem of hard-core bosons hop- 
ping on a lattice has been previously put forward by Ziegler [15] and Fresard 
[16]. Ziegler [15] has demonstrated that slave boson (SB) representation allows 
one to describe the dynamics of Bose gas taking into account three-body and 
higher effects at arbitrary densities for the three-dimensional (3D) case. The 
advantage of using slave boson representation is by the fact that there are 
only two states per site in a hard-core system: a lattice site is either empty 
or occupied by a single boson. A hopping process of a boson appears as an 
exchange of an empty site with singly occupied site. 



The paper is outlined as follows. In section 2 the Bose gas with hard-core 
interaction defined as an effective functional integral form is summarized in 
the slave boson formalism. The SB equation is evaluated in section 3. Some 
concluding remarks are presented in section 4. 
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2 The Slave Boson model 



A Bose gas, defined as a grand canonical ensemble of bosons, can be described 
in a second quantization formalism, for instance using a functional integral 
representation [18] as following: 

Z = J exp(-S[^0]/ft)nrf[0*M0] (1) 

where the Euclidean action S[<fi* , 4>] of the complex field needs to be defined 
accordingly to the model system of interest {e.g. GP, SB, ect.). In this work 
I will resort to the slave boson representation of hard-core bosons which has 
been exhaustively defined in the literatures [13,15,17,19]. The representation 
relies on the Bose-Hubbard model where the particle trades its position with 
an empty site on the lattice. Both particle as well as the empty site are de- 
scribed by the corresponding creation and annihilation operators. In the slave 
boson representation for instance, a lattice site x that is occupied by boson is 
represented by complex field b x and a site that is empty by e x . Thus the slave 
boson action S s b can be translated into the following picture [13,14] 

S s b = K e xtx,x>b x >e* x , + V x b* x b x + i\ x (e* x e x + b* x b x - 1) (2) 

x,x' 

where the first term describes the exchange of bosons and empty sites in a 
hopping process at sites x and x' with rate t XjX >. In the second term I have 
introduced an external confining potential V x . The field X x in the last term 
of the above equation enforces a constraint e* x e x + b* x b x = 1 which guarantees 
that a site is either empty or singly occupied. It is indeed the core condition 
of the slave boson model. With this choice of action, the functional integral in 
Eq. (1) can be written explicitly describing the partition function of the grand 
canonical ensembles of bosons as: 

Z = I e- s ° b l[d\ x db x db x de x de x . (3) 

The slave boson field e x and b x can be combined to form a collective field 
b* x e x — > Then the constraint field X x and the slave bosons fields can be 
integrated out which finally leads to the action for the collective field <§> x . This 
was demonstrated in detail in [15,17]. Here, I only present the results: 

Z= [ e- Sb - Sl Y[d$ x d$* x (4) 

J X 

with the hopping(" kinetic") term 

s b = Y <Mi - t)^K « £(i - r)- 1 !^! 2 + _ ?_ $ 8 (v 2 <&;) (5) 
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and the potential term 



Si = - 



2 log e 



-v x - 



V 4 j d<pexp[-tp 2 ] 



siahyf(<p- V 2 ) 2 + \$J 2 




\4 = x 2 is the planar harmonic confinement where motion of atoms are in 
planar direction with dimensionless radial distance x. The density of the con- 
densate in the trap can be calculated by the saddle-point approximation. In 
this work, we neglect the thermal fluctuations which is not important at zero 
temperatures for Bose-Einstein condensation to occurs in the system. How- 
ever, if one needed to account for the thermal fluctuations, one could do so by 
studying the deviations around the saddle point in the functional integral of 
Eq. (4) [18]. That means, we have to look for the solution order parameter <& x 
which minimizes Sb + Si taking a saddle-point approximation on the functional 
integral Eq. (4) to obtain the generalized non-linear Schrodinger equation [17] 



where a = 1— r. The term dSi/d\<S> x \ 2 contains higher-order interaction (three- 
body interactions, etc.). However this equation is difficult to solve as it is 
because the term dSi/d\& x \ 2 can lead to heavy non-linearity. 



3 Density functional solution to the slave boson equation 

Taking into account the possible effects of higher order interactions and to 
avert strong non-linearity in Eq. (7), I find a much feasible way to solve it by 
mapping to the known Kohn-Sham [20] type equation via Density functional 
theory (DFT) scheme. The main idea of the DFT formalism is in finding an 
equilibrium ground state energy by minimizing the following energy functional 



which is a unique functional of density n(x) = \& x \ 2 for a given external 
potential V x [21]. The first term in Eq. (8) is the kinetic energy of the non- 
interacting system while the second term represent the excess contribution 
of free energy E ex [n(x)}. The remainder term is the interaction energy of the 
particles in an external field. Minimizing E[n(x)], the ground state energy of 
the system can be determined by imposing the stationary equlibrium condition 







5E[n(x)} 
Sn(x) 



= V 



(9) 



4 



where the Lagrange multiplier p can be identified as the chemical potential. 
This is equivalent to the condition of equilibrium for a non-interacting particles 
under the influence of an effective external potential [22] 

8E ex [n(x)} 

v eff [n(x)] = - E ^ r + V x . (10) 

It is evident that the equilibrium density profile n(x) can be obtained for a 
given v e ff[n(x)] by solving the following Kohn-Sham type equation [20], 



5T\n(x)] r , 



$x = H$x • (11) 



By setting a = —/i and r/6a 2 = —1, Eq.(7) can be mapped to Eq. (11) 
yielding , 

p 2 S 1 [n(x)} = V H [n(x)\ + E xc [n{x)} + f dxn(x)V x . (12) 

In this relation, I have introduced an exchange-correlation energy E xc [n(x)} by 
explicitly splitting from the excess free energy E ex [n(x)} the following Hartree 
energy [22], 

Vff[n(rr)] = - J j dxdx'V(x — x)n{x)n{x) (13) 

with a generic inter-particle potential V(x — x'). In general, the exchange 
correlation functional energy E ex [n(x)} in a DFT calculation is not known 
exactly. Despite, one can resort to approximations such as the Local density 
approximation (LDA) which reads, 

E xc [n{x)\ ~ / dxE%r\p]\rMx) (14) 

where E^° m [p] = pe xc [p] and e xc [p] is the exchange correlation energy of a 
homogeneous system with uniform density p. The Functional derivatives of 
the above relation can be written as [23] 

iynM] = %^=%!*U, w . (15) 
dn{x) op 

Hence Eq. (7) can be written explicitly as the following Kohn-Sham type 
equation: 

'-V 2 + r 2 + V H [n(x)]+V 3X [n(x)]]$ x = t j,$ x , (16) 

where Vff[n(x)] = An(x)/\ \n(n(x)a 2 )\ is the Hartree potential calculated from 
Eq. (13) using a hard-core potential [7,26]. The homogeneous exchange corre- 
lation energy e xc [p] can be approximated to the following expression [25]: 

€xc [PI ~ Ml/pa*) { Hl/pa*) ) ■ {U) 
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Very recently Pilati et. al [31] did a variational Monte Carlo (vmc) calculation 
in obtaining the ground state energy of 2D bose gas using various inter-atomic 
potentials. The following is the fit of their beyond mean-field ground state 
energy : 



Fig (1) illustrate these homogeneous exchange correlation energies for vari- 
ous values of gas parameter pa 2 . Setting E xc [n(x)] = one recovers the GP 
equation. 

There are three cases to consider in solving Eq.(16), namely using information 
e x° l °lp\i e imc[p\ or E xc [n(x)\ = 0. The KS equation Eq.(16) need to be solved 
self-consistently with the condition that the areal integral of the condensate 
density n(x) = \& x \ 2 is equal to the total number N of particles. For a numer- 
ical illustration, we have taken the radial trap frequency, and the scattering 
length as appropriate for 23 Na atoms in the experiment of Gorlitz et al. [1], 
namely u>± = 188.4 Hz, s-wave scattering length a = 2.8 nm. I am implic- 
itly assuming however, that the trap has been axially squeezed to reach the 
strictly 2D scattering regime. In this work all the distances and energies have 



by fkv±/2. Numerical iteratation method by discretization using a two-step 
Crank-Nicholson has been used in performing this calculation for 5 x 10 3 and 
5 x 10 4 23 sodium atoms. 

The results obtained from Eq. (16) using information of the exchange cor- 
relation energies e x ° l [p] [rikoi(x)} and e^™ c [p] [n 1)mc (a;)] comparing to the case 
with E ex [n(x)] = [^gp( 2; )] are depicted in Fig. (2) and Fig. (3) for 5 x 10 3 
and 5 x 10 4 atoms respectively. The main effect of the correlation potentials 
can be observed where the condensate density reaches the peak at the trap 
center. So it is feasible to explain the observation by measuring the differences 
An* = nfcoi(O) — ^gp(O) and An = rikoi(Q) —n vmc (0) at the trap centre as done 
by Nunes [24]. For a fixed number of atoms N the conditions / dxAn* = 
and / dxAn = are imposed so that the rest of curves An*, An are well de- 
termined at x = 0. The differences An*, An increases with An'/nfoi, An/n^ 
varying from 9% and 1.33% for TV = 5 x 10 3 reaching about 13% and 2.27% for 
N = 5 x 10 4 respectively. The increment in An*, An corresponds to the mag- 
nitude of reduction in the effective interaction energy of the system. The small 
increment in An shows that the density profile predicted using the correlation 
energies e v ™ c agrees well with the one obtained using e x ° l . However what is 
more spectacular is the large values of An* and the increment of An/n^ from 
9% to 13% as the number the number of atom increases from N = 5 x 10 3 to 
N = 5 x 10 4 . The effect of exchange-correlation clearly shows a reduction in 
effective interaction energy thus the profiles are more localized at the centre 
of the trap in comparison to the GP case. This argument can also be sup- 
ported by the observation in Fig (1) where the magnitude of homogeneous 




(18) 



been scaled in the harmonic oscillator 
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exchange correlation energies decreases as the density increases. On the other 
hand e xc [p] is negligible for the very dilute regimes pa 2 <C 0.01 justifying GP 
approximation is only a good model in this limit 



4 Conclusion 

In summary, the harmonically trapped 2D Bose gas at zero temperature de- 
scribed by the slave boson model has been studied for both the dilute and 
dense limits. The central issue of this work is in how to treat the slave boson 
term dSi/d\<$> x \ 2 in the highly non- linear partial differential equation (7). The 
importance of the higher order interactions via the exchange-correlation en- 
ergy has been revealed by mapping Eq. (7) to the Kohn-Sham scheme of DFT 
in the dense limit. This latter term is responsible for the correction of up to 
9% and 13% of the peak density for iV = 5 x 10 3 and iV = 5 x 10 4 atoms 
at the center of trap. In other words these exchange correlation has restored 
the correct high density behaviour of the inhomogeneous Bose condensate im- 
posed by the strict hard-core condition of the slave-boson representation. On 
the opposite dilute regime I have demonstrated that the slave boson equation 
(7) can be approximated to a Gross-Pitaevskii type equation (via DFT). 

At absolute zero temperature where the thermal fluctuation is frozen out and 
quantum fluctuation is prominent. This quantum fluctuation plays an impor- 
tant role in the phase transition from superfluid to mott insulator [27-29] in 
the three-dimensional case. Recently Wu et.al. [30] have demonstrated that 
a vortex in 2D optical lattice can induce the SF-MI transition. As a future 
direction a practical realization of this transition will be investigated using 
the slave-boson model of the 2D condensate accomodating a single or many 
vortices. 
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Fig. 2. Density profiles n(x) for the condensate (in units a^N- 1 , with 
dho = \J h/muj±) versus the radial distance x (in units of a^ c ) obtained by solving 
Eq.(16) using information of e^°' (solid line), e v ™ c (long-dashed lines) and setting 
E xc [n(x)] = (short-dashed lines). 
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